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Bistatic Scattering From Conducting 
Calibration Spheres 


1. INTRODUCTION 


After an extensive literature search, no comprehensive handbook of scattering 
tables for conducting spheres was found, Since a handbook of this kind would be 
valuable for calibration procedures, such as the use of spheres for the scattering 
reference targets in anechoic chambers, a project was undertaken to develop this 
handbook. 

Mie! has generally been recognized for deriving the basic formulas to calculate 
the complete analytical solution for the scattering from spheres of arbitrary size. 


However, a literature review conducted by Logan shows that the formulas con- 
tained in Mie's solution were first presented by Lorenz* in 1890. One reason why 
Mie has been credited with the solution may be that he set out on an extensive pro- 
gram of applying "his" solution to the computation of scattering from spheres which 
were too large to fit the criteria for Rayleigh scattering. In spite of the fact that 
Mie's solution was not original, it is exact and is used as the basis for the computa-~ 


tion of the scattering tables contained in this handbook. 


(Received for publication 18 April 1984) 
1, Mie, G. (1908) Ann, d. Physik, ae -442, 


2. Logan, N.A. (1962) Early history of the Mie solution, J. Opt. Soc. Am., 
52:342-434, 


3. Logan, N.A, (1965) Survey of some early studies of the scattering of the plane 
waves by a sphere, Proc. IEEE, pe a 


4. Lorenz, L. (1890) Videnskab. Selskabs Skrifter, aK 


. 


ee a 
ny z 
a : 
a A number of existing FORTRAN programs for calculating the Mie scattering for oe 
* spheres were analyzed for completeness, accuracy of computation, and ease of Pe 
, adaptation. 5, 6.7 We finally chose the program developed by Hancock and k 
- Livingston’ of the Naval Research Laboratory (NRL) because of: (1) its complete an ve 
- and accurate treatment of the Bessel function computations as described in Abramo- my 
* witz and Stegun;® (2) its computational assurance that each Mie series reaches a os 
designated convergence; and (3) its convenient adaptability for use on the Air Force : hoy 
a7 Geophysics Laboratory (CDC 6600 computer). , 
ee Each table in this handbook lists the radar cross section normalized to the we 
_ geometric cross section for both E-Plane (E)) and H-Plane (E.) polarizations for oe 
eo bistatic scattering angles from backscatter to forward scatter at 2-degree intervals. a 


A table is listed for each ka value ranging from 0.2 to 40 at intervals of 0.2, = 
In addition to the tables, two sets of plots are included in this handbook. These 
plots are intended to illustrate trends and to give the observer an overall perspective pane 


S of the bistatic scattering from conducting spheres. The first set contains plots of a 
. Ofn a2 vs ka for particular bistatic observation angles ranging from backscatter to ae 
forward scatter at 10 degree intervals. The second set contains plots of ofnae vs aa 
“4 bistatic observation angles for selected ka values ranging from 0.01 to 70. Ly 
S 2. THEORETICAL FORMULATION as 
7 The FORTRAN program is written from Kerker's formulation? of the Mie scatter- | 
- ing solution. A flowchart is contained in Appendix A and the program itself is listed oe 
i in Appendix B. The NRL version of this program computed the scattering from loss- ix 
less dielectric spheres only and, therefore, required modification before computa-~ he 
es tions for perfectly conducting spheres could be carried out. The derivation of the Nn 
formulas for conducting spheres will now be discussed. i 
@ The geometry of the problem is illustrated in Figure 1, where a sphere of om 
rs radius a is centered at the origin. The magnitude of the incident electric field is oa 
as assumed to be one propagating in the +z direction, én and polarized along the x-axis ap 
a (E;). The observation direction of the scattered wave (r) is defined by the usual 7 
spherical angles @ and @ shown in Figure 1. The directions of incidence and scatter - ; 
- ing define the scattering plane." : 
4 = 
e) mS 
: (Because of the large number of references cited above, they will not be listed here, ed 
on See References, page 55, ) : oy 
a = 
; 7 f 


a, as 


The most difficult part of calculating the radar cross section (RCS) is the 
evaluation of the complex scattering amplitude functions S 1 and So given by 


J 
es 2n+1 
s, = > nin 1) [ ah 7 (cos 6) + db, a (cos 6) | (1) 
2 2n+1 
So = Leal aot [ a, T,,(cos 6) + db. 7 (cos 9) » (2) 


where the angular functions 7, (cos 6) and t “ (cos @) are 


pil) (cos 9) 
n 


7 (cos @) = ain (3) 
T , (cos @) = a BO cos 6. (4) 
From Eq. (3) 
p (cos 6) = sin @ * 7 (cos 6) (5) 
which leads to 
T , (cos 9) = 5 p! (cos 6) = cos @ * m, (cos @) - sin? 9 ° 7 (cos 6). 
(6) 


The following recursion relations can be used to compute 7, (cos 6) and 7 (cos @). 


7, (cos 6) = cos 6 (223) Tm, - 1 (eos 9) - (2) T - 2's 9) (7) 
m. (cos 8) = (2n-1) m,_, (cos @) + m!_ 9 (cos @) (8) 


knowing that 7, (cos 6) =1, 7 (cos 6) = 0, 1, (cos 6)= 3 cos @, and 1, (cos 6) = 3. 


From the scattering amplitude coefficients, an and bw’ for penetrable spheres, 


ET EE Te et EO ee Ee ee Le ee ee eee OY 


3”. 


) 


7? OST 


rs 


wee 
° 
e 


(9) 


e ‘ 
= Hn (0) Hn 8) - mH, (8) ¥, (ad 
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we can find an and b, for a perfectly conducting sphere by letting the refractive 
index (m) approach infinity. Thus Eqs. (9) and (10) reduce to 


a Ser (11) 
mn & (a) 
and 
ve (a) 
bo = tT . (12) 


The perimeter of the sphere in wavelengths is denoted by @ = ka = 203. with r= free 


space wavelength of the incident wave and a = radius of the sphere. ¥ (a) 
and Cm are the Ricatti- Bessel functions defined in the following manner: 


1/2 
WO = (ala J ay (@ (13) 
Gta) = (nasal? yl?) ay, (4) 


(2) 
n+1/2 
first kind and Hankel functions of the second kind respectively. The function, 


where Jne 1/2) and H (a) are the half integral order Bessel functions of the 


¢,(@) is computed fium ¥ (a )~ ix ,(@) where 


x,(@) = (na/2ytl2 y (15) 


n+1/2 (a) 
and aw 1/2') is the half integral order Bessel function of the second kind. Primes 
on the functions ve (a) and ie (a) represent the derivatives with respect to the 


argument (a). 
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ou 
The method of logarithmic derivatives, first used by Infeia!®, will now be oe 
used to transform Eq. (11) into a computationally simpler form which does not ws 
involve the explicit evaluation of vita) and cr (a). We now define two Bessel ‘. ot 
density functions S,(@) and R(@) where a 
a 
1/2 hoe 
7 1/2 _ [2 eS 
S,(@) = (a) Jeg 1/2!) = (3 ) via) (16) : — 
1/2 Hie 
5 1/2 ,,(2) _ {2 ee 
Ri(@) = (a) H+ 1/2!) = (2) Sn (a). (17) He 
The logarithmic derivatives of Vv (@) and fa) are defined as two new functions oa 
se: 

(a) and p fa): - 
a ¥i(a) Sita) sg 

a (a) = aa In Wo) = ¥_tay = Stay~ Ja in S,'@) (18) 
a tia) Ria) g = 
p. (a) =-~ In | & (a)} = = = In | R (a) . (19) 
n da n Tay Ray da ne |) nS 
The expression for the scattering amplitude coefficient a, now becomes oe 
, — 
' ¥ (a) a i 
wa) Ya" Fe via) oa) Pas 
“nar © tian aT ae 
ola) + ee 
< cfm) a 
By substituting Eq. (12) into Eq. (20) we get nee 
g(a) on 
a= b, . pay ae (21) oo 
If we let z (a) denote the general Bessel function [J (a), ¥ (@), H(q@), and so on], — 
then the general Bessel density function D(a) becomes yo 
= gil? aoe 
Dia)=a Zane 1/2'%) > (22) oe 
ee 
che 
10. Infeld, L. (1947) The influence of the width of the gap upon the theory of ee 
antennas, Quart. Appl. Math., 5. . sae 
eter ee 
Sea 
1 iu 
= 


By differentiating D(a) with respect to a we get 


e a 1/2 ‘ 1 -1/2 
11 
If the recurrence formula 
Zao 1/26 = 2-1/2! a Z 41/2! (24) r 
is substituted into Eq. (23) then 
Pie Ae _n Bae 
Dia) = a & 12 - B24 2| (25) 


If we let d (a) denote the logarithmic derivative of Da) then 


D! (a) Zz (a) 
.~ _n = n~1/2 
d(a) = dD,” 


-2 (26) 
n+ 1/2 o 


and the computation of O(a) and Pp (a) become 


wi (a)  S* (a) J (a) 
ota) = Py = stay = peor _ 2 (27) 
n ni n n+1/2 iy 
’ £ (2) 
ee C(@) Bylo) . He 1/2 (@) a ae 
n tT Ta) R, (a H a * 
z n+1/2(a) 
By substituting Eqs. (13) and (14) in Eq. (12) we find 
¥,, (a) nar/2yi!® Ing 3/2) - : +1/2'%) 
== (29) 


Be oa le ee ay ; 
e n'®) (pa /2)'° H (a) H 9 (a) 


s : 
te 


n+1/ 


a 
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11. Watson, G. (1962) A Treatise on the Theory of Bessel Functions, Cambridge 
University Press, New York. 
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These relations are sufficient to easily compute the RCS of perfectly conducting 
spheres and were utilized for modifying NRL's FORTRAN program. Because we 
preferred the representation and notation of Ruck et al, 12 we decided to convert 
Kerker's formulation to Ruck's notation. Figure 2 shows the geometrical relation- 
ship between the two different notational systems. As one can see, Ruck's o” is 
equal to Kerker's 1/2 - @ and Ruck’'s @ ’ is equal to Kerker's 7 - §. After this 
conversion the normalized RCS, one which is the absolute RCS in the observation 
direction divided by the geometric cross section of the sphere, becomes 


2 
oy (0', 0’) = (=) [!s, (9')|? cos” 9’ + |, (0")|? sin” 0"| (30) 


3. PROGRAM DESCRIPTION 


The first modification of NRL's computer program was the conversion to 
CDC 6600 single precision arithmetic (approximately 14 decimal digits). Compari- 
son of single and double precision runs showed no difference for the four decimal 
digit accuracy desired. Therefore, double precision was considered unnecessary. 
As in the original program, the convergence criterion was that six consecutive terms 
of each series are required to satisfy the condition that the ratio of the current 
term to the current partial sum not exceed the relative tolerance of + 0.001. If this 
convergence condition is not met within 90 terms, a convergence warning is issued, 
and the 90th partial sum is used as the series value. Warnings have only been ob- 
served for ka > 77. The NRL program was further modified to be more user 
friendly. In other words when a data input is required an explicit question appears 
as aprompt. The NRL program only allowed the user to compute the RCS at 19 
equispaced bistatic scattering angles but our modified program allows an unlimited 
number of equispacei angles to be computed in a single run. Subroutines in the 
NRL program which were used solely to compute the scattering from particle dis- 
tributions have been omitted. Other changes include the conversion to Ruck's 
notation, the change in the calculation of scattering amplitude coefficients an and b, 
to accommodate perfectly conducting spheres, and expansion in capability to compute 
the normalized RCS, aN for any two orthogonal polarizations (not just E-Plane 
and H-Plane). 


12. Ruck, G., Barrick, D., Stuart, W., and Krichbaum, C. (1970) Radar Cross 
Section Handbook, 1, Plenum Press, New York. 
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*, There is a choice of two input formats. Type 0 format allows for ka and A as 

input constants; type 1 allows for frequency and sphere radius, For either type, 

the normalized RCS, 9,» for two orthogonal polarizations @’ and ¢’ + 90 is cal- 

2 culated for scattering angles 9’, 9'+ 9", 9'+ 240’, ... 6’ + (NRTH - 1)40’. : 
The input variables are defined below: 


LAMBDA A Free space wavelength 
A a Sphere radius (centimeters) 


27a 
KA ka = 


FREQ f Frequency (Hertz) 
PHI ¢ Labeled PHI’ - Polarization angle in XY plane (degrees) 


THZ Labeled THETA’ SUB ZERO - Initial value of @’ (degrees) 


“ DELTH 46 Labeled DELTA THETHA’ - Step size for 6’ (degrees) 


NRTH Labeled NUMBER OF THETA’ - Number of values of 6% 
to be computed 
The first line of output for either type 0 or 1 run contains NRTH, A (meters), 
LAMBDA (meters), and ka respectively. Each successive line contains 6’, 
: Oy (8 Vs Oy (8 ’, 9’ + 90°). Two sample runs are listed on the next page. 
The first is for type 0 runs; the second, type 1. 


4. PLOTS OF 0/Ta2 VERSUS ka 


This section contains plots of o/na” vs ka for conducting spheres. Pilots for 
each bistatic scattering observation angle ranging from backscatter (@ = 0°, 
where @ corresponds to Ruck's 6) to forward scatter (9 = 180°) at 10 degree inter- 


ae vals are included. The solid curves represent the E-Plane scattering and the 


re dashed curves represent the H-Plane scattering. When referring to these plots, 


. 
' 


three distinct regions will be identified: the low-frequency region where ka < 0.4, we 


ee the resonance region between ka = 0.4 and a point where the oscillations have 


? 


settled to within some smal! value, ¢€, of the asymptotic value of o/nae and the 
high frequency region where o/naz is within € of the asymptotic value. woe! 
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INPUT DELTA THETA , NUMBER OF THETA’ *10,19 
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TYPE Q FOR KA, 1 FOR FREQ INPUT >I 
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6.1788E-01 
3.7935E-91 
2.8147E-01 
3.3200E-01 
5.1625E-01 
? 9635E-01 
1. 1157E+00 
1,4087E +00 
1.61386 +00 


180.0 1.6875&+00 


3.6116E+90 
3.5760F +00 
3.5209E +00 
3.4425E+00 
3. 5376E +00 
3.2046E +00 
3.0448E +00 
2.8628E +00 
2.6664E+00 
2.4661E+00 
2.2732E+00 
2.0986E +00 
1.9506E+00 
1.8345E+00 
1.7522E€+00 
1.7035E+00 
1.68756 +00 


20.0 
30.0 
40.0 
50.0 
60.0 
70.0 
80.0 
90.0 
100.0 
110.0 
120.0 
130.0 
140.0 
150.0 
160.0 
170.0 
180.0 


3.4103E+00 
3. 1361E 400 
2.7732E +00 
2.3434E +00 
1.8759E+00 
1.4058E +00 
9.7337E-01 
6.1948E-01 
3.B078E-01 
2.8279E-01 
3. 3331E-01 
5.9771E-01 
7, 9811E-01 
1. 1178E+00 
1.4112E+00 
1.6166E+00 
1.6904E+00 


-2998E-01 
3.6385E+00 
3.6322E+00 
3.6127E+00 
3.5775E +00 
3.5225E +00 
3. 4443E+00 
3.3397E+00 
3.2069E +00 
3.0474E+00 
2 .8655E +00 
2.6693E+00 
2.4690E +00 
2.2761E+00 
2.1015E+00 
1. 9535E +00 
1.8374E+00 
1, 7551E+00 
1.7065E +00 
1,.46904E+00 


1.001 


SOREEEESRORESE EHEC ETOH 
TYPE 1 TO RUN, O TO STOP v6 
ENR MIECH 
035500 MAXIMUM EXECUTION FL. 
.097 CF SECONDS EXECUTION TIME. 


DRECHEEEEEESESETESERSRE 
TYPE 1 TO RUN, O TO STOP 20 
END NIECH 
035500 MAXIMUM EXECUTION FL. 
.099 CP SECONDS EXECUTION TIME. 
COMMAND- 


COMMAND- 
Figure 3. Sample Run for MIECM 
for Ka, Lambda Inputs 


Figure 4, 


Sample Ran of MIECM 


for Frequency, Radius Inputs 
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Figure 5 shows the classic normalized backscatter cross section vs ka which 


is seen in most reference books on electromagnetic scattering. The remainder of 


P the plots included in this section are rarely seen in the literature, but they do ; 
: reveal some trends which should be noted. First, as @ increases from 0 degrees, ‘ 
*. the resonance region for the E-Pliane scattering broadens as the oscillations 

* increase in magnitude. Thus, as @ increases, larger ka values must be attained S 


before reaching the high-frequency region. On the contrary, for @ < 80°, the: 
en H-Plane scattering curves remain relatively unchanged and reach the high-frequency 
, region at much lower ka values than for the E-Plane. Second, the H-Plane curves 
indicate a forward scatter enhancement as @ increases beyond approximately 
: 80 degrees. This enhancement does not become apparent for E-Plane scattering 
until 6 reaches much higher values. Third, the E-Plane and H-Plane scattering 


_ curves clearly begin to reconverge as 6 increases beyond 140 degrees. Fourth, 
at low values of both 6 (except 6 = 0°) and ka, the E-Plane scattering curves 
ay indicate the presence of a second order interference pattern. This resonance 


manifests itself as a modal pattern in those E-Plane curves. The H-Plane curves 
reveal no such phenomenon. Lastly, for the forward scatter direction (Figure 23), 
it is clear that for ka > 10 we find that o/na2 = (kay? which agrees with the high- 
frequency region solutions for forward scatter universally found in the literature. 


5. PLOTS OF o/ a2 VERSUS 6 


Si This section contains plots of o/nae vs the bistatic scattering observation angle 
6 (where @ corresponds to Ruck's 6 ‘) for selected values of ka ranging from 0.01 to 
70,0, The solid curves are for the E-Plane scattering and the dashed curves are 
for the H-Plane scattering. 

Figures 24 and 25 illustrate the dipole effect which is prevalent in very small 
On conducting spheres. The induced surface charges and currents on such spheres 
: correspond to oscillating electric and magnetic dipoles where the ratio of the elec-~ 
& tric to magnetic dipole moment is 2 to 1. The null in the E-Plane curves at @ = 120° 
; is a result of the electric and magnetic dipole far fields cancelling at this angle. 
= Specifically, the low-frequency scattering solutions are 
ss 2 


o/na 4 (ka)? {1/2 + cos 6}° for E-Plane 


= g/na? 


4 (kay* [1/2 cos 6 + ‘ie for H-Plane 


which clearly reveal the null in the E-Plane. It is also evident from the first two 
plots and the formulas above that o/na” for H-Plane scattering is a monotonically 


decreasing function containing no nulls. 
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As ka increases beyond 0. 2, the null in the E-Plane curves gradually dis- 
appears and no clear trend appears until ka increases beyond approximately 4.0. 
The envelope for the E-Plane curves then becomes evident and even though the 
broadness and magnitude of the envelope pattern varies as ka increases, the 
characteristic shape does not. For the H-Plane, the magnitude of the oscillations 
is always much smaller. 

One important factor to consider when using metal spheres for calibrating 
bistatic scattering ranges is that one must be careful not to pick a sphere size 
which has a deep valley or null in the theoretical pattern at the bistatic angles to 
be measured. This condition must be insured in order to minimize measurement 


errors throughout the dynamic range of the measurement system. 


6. CONCLUSIONS 


We have illustrated the theoretical derivation for the Mie scattering from con- 
ducting spheres using Kerker's formulation and have incorporated that theory into 
the FORTRAN program listed in Appendix B. This program (MIECM) is capable 
of computing the normalized scattering cross section for any number of equi- 
spaced bistatic scattering observation angles for any two orthogonal polarizations 
(for example, E-Plane and H-Plane). A description of how MIECM operates is 
provided so that it can be used after a minimal amount of user familiarization. 

Two sets of plots which illustrate trends have been included in this handbook, 
In the discussion some anomalies such as the presence of a second order inter- 
ference pattern in the E-Plane scattering for low values of @ and ka, and the dipole 
effect which is prevalent in the E-Plane scattering of very small spheres are pointed 
out. It is hoped that these two sets of plots will give the observer an overall 
perspective of bistatic scattering from conducting spheres. 

In addition, we have stressed that, in order to minimize measurement errors, 
the size of the sphere used in the calibration of a bistatic scattering range must be 
chosen very carefully. 

Finally we have compiled and included a comprehensive set of bistatic scatter- 
ing tables for perfectly conducting spheres. The tables are intended to serve as a 
reference for the calibration of indoor scattering ranges which employ conducting 
spheres as test bodies. These tables will minimize the need for computing the MIE 


scattering for most reference spheres, 
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Appendix A 


Flow Chart of Program MIECM 


PROGRAM MIECM 


IF ITYPE-0: STOP 


READ DATA 


COMPUTE ANGLE FUNCTIONS 
7 ,(cos8) AND 1, (cos@) FOR ALL 6 


CALL DSPHJN AND DSPHYN 
COMPUTE TABLES OF SPHERICAL 
BESSEL FUNCTIONS j,, AND y,, 


FORO <n<90 


CALL ANBN (IONETWO) 
IONETWO = 1: 
COMPUTE a, AND b,, 


CALL SCAT 
COMPUTE SCATTERING AMPLITUDE 
FUNCTIONS |S,|2 AND!S,|2 

FOR ALLO 


COMPUTE ADDITIONAL a, AND by 
AS AND IF NEEDED FOR S, AND So 


CALL STEP 
STEP THE 4 SERIES FOR S, AND S, 


CONVERT TO RUCK’S NOTATION 


PRINT RESULTS 
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Appendix B 


phy Program MIECM Listing 


PROGRAM MIECM( TAPE1, TAPE2, TAPEG) 000100 


” COMMON/ALL/FACTR.TOL,NTOL,KOUT, PN(19,90),TN( 19,90) 000110 
. COMMON/SCAN/ALF , ANI (162), ANR(162).8Ni(162),8NR(162), 000120 
i. 1 CSC,CXT, LMDA,NDEL,NMAX,PI,PNPR( 19,3) ,NCSC,NCXT 000130 
7 TYPE REAL KA,LNDA 000140 
ey COMMON/MSC/A,DELTH, ITYPE,NMAXSAV,NTOP, THZ 000150 
we COMMON/MS I SC/CAPKSC , CAPKXT ,$1S0(19) ,S2SQ(19) 000160 
= DATA( TEST=0.005),(LIM#t0) 000170 
cs REWIND 2 000180 
in CALL CONNEC(1) 000190 
CALL CONNEC(6) 000200 

P1=3.1415926535898 000210 

s NMAX SAV2NMAX 000220 
= PRINT(6, 18) 000230 
a 18 FORMAT (24(1H®*)) 000240 
= c 000250 
a C INPUT INITIALIZING INFO 000260 
. c 000270 
‘ 20 PRINT(6,*)“TYPE 1 TO RUN, O TO STOP >* 000280 
a READ (1,#) ITYPE 000290 
IF(ITYPE .£Q. 0) GO TO 67 000300 

PRINT(6,*)"TYPE O FOR KA, 1 FOR FREQ INPUT >* 000310 

READ (1,#) LTYPE 000320 

. IF(ITYPE.EQ.0) GO TO 22 000330 
= PRINT(G.*)"INPUT FREQ, RADIUS(CM) >* 000340 
om READ (1,%) FREQ,A 000350 
LMDA=2.998E8/ FREQ 000360 

¥ A=A*.01 000370 
" KA=2.0¢PI*A/LMDA 000380 
4 GO TO 25 000390 
22 PRINT(S,*) "INPUT KA,LAMBOA(CM) >* 000400 

READ (1,*) KA, LMDA 000410 

ae LMDA=LMOA*.01 000420 
A=KA*LMOA/2.0/PI 000430 

25 PRINT(6,*) "INPUT PHI', THETA' SUB ZERO >* 000440 

READ (1,*) PHI,THZR 000450 

PRINT(6,*#) “INPUT DELTA THETA’, NUMBER OF THETA’ >* 000460 

READ (1,*) DELTHR,NRTH 000470 

PHI1=PHI+90.0 000480 

c 000490 

C CONVERT POLARIZATION ANGLES TO RACIANS AND FROM KERKER'S 000500 

C COORDINATE SYSTEM TO RUCK'S COORDINATE SYSTEM 000510 

c 000520 

e PHI=(90.9-PHI) *P1/180.0 090530 
a PHI12(90.0-SHI1)*#PI/180.0 000540 
.* CSPHISQ=COS( PHI) *COS( PHI) 000550 
3 SNPHISQ=SIN(PHI)*SIN( PHI) 000560 
3 CSPH1ISQ2CNS(PHI1) *COS(PHI1) 000570 
’ SNPH1SQ=SIN(PHi1)*SIN(PHI1) 000580 
’ LPHI=TFIA( PHI) 000590 
IPHIVSIFIX(PHI1) 000600 

c 000610 

C CONVERT SCATTERING ANGLES FROM RUCK'S THETA' NOTATION 000620 

. C TO KERKER'S THETA NCTATION 000630 
. C SEE REFS. “RADAR CROSS SECTION HANOBOOK, VOL 1" BY RUCK 000640 
. C ,ET AL, AND “THE SCATTERING OF LIGHT AND OTHER 000650 
: C ELECTROMAGNETIC RADIATION* BY KERKER FOR DETAILS 000660 
. c 000670 
. THZ2180.0-THZR 000680 
; DELTH=-DELTHR 000690 
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24 
26 


c 


PRINT(2,26) NRTH,A,LMDA,KA 
FORMAT(14,€11.4,£11.4,F8.3) 
FNRTH=FLOAT(NRTH) 
PM1=2—NRTH/19- 
MP1=INT(PM1) 
FMP1=FLOAT (MP1) 
MNRTH=MP1—-1 
IF(PM1 .GT. FMP1) MNRTH=sMP1 
NNRTH=NRTH-19°MNRTH 
MNRTH=MNRTH+ 1 
DO 66 MRTH=1,MNRTH 
NTOP=19 
IF(MRTH .€Q. MNRTH) NTOPsNNRTH 
00 26 N=1,NTOP 
DTS=(MRTH-1) *194+N=1 
TH RAD=(TH Z+0TS*DEL TH) +P1/180.0 
COS TH=COS(TH RAD) 
SIN TH=SIN(TH RAD) 
PN(N,1)=1.0 
PN(N,2)23.0*COS TH 
TN(N,1)=COS TH 
TN(N,2)23.0°(COS TH**2=SIN TH#*2) 
PNPR(N,1)=0. 
PNPR(N,2)=3.0 
DO 24 Kz3,NMAX SAV 
OKzK 
TS=2.0*DK-1.0 
Run mieces Fue (OTS/ (DMS 1.20) 1ePN(N KA) >(DK/ COWL: O))SENCM SRE 2) 
PNER(N, 3) eD0TS*PN(N,K=1)+PNPR(N, 
TN(NeK)2COS THePN(N,K)“SIN TH*SIN THePNPR(N,3) 
PNPR(N,1)sPNPR(N,2) 
PNPR(N, 2) =PNPR(N, 3) 
CONTINUE 
CONTINUE 
CALL SCAT 
00 36 I=1,NTOP 


C CONVERT SCATTERING ANGLES BACK TO RUCK'S THETA' NOTATION 


c 


34 
36 
66 


67 


aan 


THSNGL=THZR+( (MRTH=1)*19+1-1) sDELTHR 
TS=(LMDA/PI) #82 

PANT terse (51SQ( 1) CSPHISQ+S250( 1} *SNPHISQ) 
PRNTI=PRNT1/As#2 
PRNT2=TS*(51SQ( 1) *CSPH1SQ+S25Q( 1) *SNPH1SQ) 
PRNT2=PRNT2/As#2 

PRINT(2,34) TH SNGL,PRNT1,PRNT2 
FORMAT(F6.1,1P,2E11.4) 

CONTINUE 

CONTINUE 

PRINT(6,18) 

GO TO 20 

CONTINUE 

CALL DISCON(1) 

CALL DISCON(6) 

REWIND 2 

END 


SUBROUTINE STEP(I,K,PN FACT,TN FACT ,SUM, TERM,NCONSEC,N USED) 
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000700 
000710 
000720 
000730 
000740 
000750 
000760 
000770 
000780 
000790 
000800 
000810 
000820 
000830 
000840 
000850 
000860 
000870 
000880 
000890 
000900 
000910 
000920 
000930 
000940 
000950 
000960 
000970 
000980 
000990 
001000 
001010 
001020 
001030 
001040 
001050 
001060 
001070 
001080 
001090 
001100 
001110 
001120 
001130 
001140 
001150 
001160 
001170 
001180 
001190 
001200 
001210 
001220 
001230 
001240 
001250 
001260 
001270 
001280 
001290 
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20 
22 


aan 


26 


30 


32 
33 


waa ats! ta 
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COMMON/ALL/FACTR, TOL,NTOL,.KOUT, PN( 19,90), 1N(19,90) 
TF(NCONSEC .E€Q. NTOL) RETURN 

TERM=FACTR*(PN FACT*PN(I,K)4TN FACTOTN(1,K)) 
SUM=SUM+TERM 

N USED=N USED+1 

1F(ABS(TERM/SUM) .GT. TOL) 20,22 

N CONSEC20 

RETURN 

N CONSECsN CONSEC+1 

IF(N CONSEC .EQ. NTOL) KOUT#KOUT+1 

RETURN 

END 


SUBROUTINE SCAT 
COMMON/ALL/FACTR, TOL,NTOL,KOUT.PN(19,90) ,TN(19,90) 
COMMON/SCAN/ALF, ANI (162), ANR(162),BNI(162),BNR(162), 
1 CSC,CXT, LMDA,NDEL ,NMAX, PI, PNPR( 19,3} ,NCSC,NCAT 

TYPE REAL LMOA 

COMMON/MS 1SC/CAPKSC , CAPKXT ,$159(19) ,S2S9(19) 
COMMON/MSC/A,DELTH, ITYPE, NMAXSAV,NTOP, THZ 

DATA( TOL#1E-3). (NTOLS6) , (NMAX#90), (NDEL#10) 

NMAX=NMAX SAV 

ALF22.0¢PIeA/LMDA 

1 ONE TwO#! 

NLO#1 

NHI =NMAX 

CALL AN BN(I ONE TWO,NLO,NHI,NTAB) 

CAP K SCaCSC/(PI*Aea) 

CAP K XTsCXT/(PISA*A) 

DO 36 121,NTOP 

SiIR=S1[=S2R=' 2120. 

N SIREN S1I3N S2R2N $2180 

N S1R TOLsN StI TOL=N S2R TOLeN S2I TOLed 

K DUT=0 

KLO=1 

KHI=NTAB 

pO 30 K=KLO,KHI 

OK=K 

FACTR=(2.0*DK+1.0)/(DK*(DK+1.0)) 

CALC STEP(T,K,ANR(K), BNRIK), S512, TERM SIR,N SIR TOL,N SIR) 
CALL STEP(I,K,ANI(K),BNI(K),S12,TERM S1I,N StI TOL,N S11) 
CALL STEP(1,K,BNR(K),ANR(K),S2R,TERM S2R,N S2R TOL,N S2R) 
CALL STEP(I,K,BNI(K),ANI(K),S21,TERM S21,N S2t TOL,N S2I) 
IF(K OUT .€Q. 4) GO TO 32 

CONTINUE 

NLO=KLO=NTAB+1 

NHI =KHI=MINO(NTAB+NDEL ,NMAX) 

IF(NHI LT. NLO) GO TO 32 

1 ONE TwO=2 

CALL AN BN(I ONE TWO,NLO,NHI,NTAB) 

GO TO 26 

MINTOL=MINO(NSTRTOL ,NSUITTOL,NS2RTOL ,NS2ITOL) 

IF(MINTOL LT. NTOL) PRINT(6,33)MINTOL 

FORMAT(* CONVERGENCE WARNING--ONLY*13¢ TERMS SMALL®*) 
S1SQ( 1) #StRee2+Stlee2 

$25Q(1) 2S2Re*2¢S2le02 

CONTINUE 

RETURN 


$4 


001300 
001310 
001320 
001330 
001340 
001350 
001360 
001370 
001380 
001390 
001400 
001416 
001420 
001430 
001440 
001450 
001460 
001470 
001480 
001490 
001500 
001510 
001520 
001530 
001540 
001550 
001560 
001570 
001580 
001590 
001600 
001610 
001620 


001630. 


001640 
001650 
001660 
001670 
001680 
001690 
001700 
0017140 
001720 
001730 
001740 
001750 
001760 
001770 
001780 
001790 
001800 
091810 
001820 
001830 
001840 
001850 
001860 
001870 
001880 
0018690 
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001900 


c 001910 
c 001920 
c 001930 
SUBROUTINE ANBN(IONETWO,NLO,NHI ,NTAB) 001940 
COMMON/ALL/FACTR,TOL,NTOL.KOUT.PN(19,90),7N(19,90) 001950 
COMMON/SCAN/ALF , ANI (162), ANR( 162), BNI (162),8NR(162), 001960 

4 CSC.cxT, LMDA, NOEL ,NMAX, PI, PNPR( 19.3) ,.NCSC,NCXAT 001970 

TYPE REAL UN. L¥DA 001980 
OIMENSION UN{ 300), YN( 162) 001990 

TYPE COMPLEX AN,6N,HN,HN1,RHON 002000 
OATA(CUTetEtzo} 002010 
IF(IONETWO.EQ.2) GO TO 18 002020 

12 CSCsCxT=0.0 002030 
NCSC2NCXATZNSCTOL=NXTTOL&O 002040 

16 CALL SPHUN(ALF,NHI, UN) 002050 
CALL SPHYN(ALF NHI, YN) 002060 

T=0 002070 

161 Ise! 002080 
1F(ABS(UN(14+1)).GT.1E5) JN( 141) UN(1I+1)/1E270 002090 

IFC ABS(YN(141)). LT. 1E-5) YN( 341) 8YN( 141) 014E270 002100 
1F(ABS( UN(14+1)).GT.1.0/CUT.AND.ABS(YN(341)).LT.CuT) GO TO 17 002110 
NMAX=NHI «1-2 002120 

GO To 18 002130 

17 IF(I1.LT.NHI) GO TO 161 002140 
18 00 44 NsNLO,NHI 002150 
ON=N 002160 
HN*CMPLX(UN( N41) ,-¥N(N44)) 002170 
BN=JN(N+1)/HN 002180 
SIGMANs JN(N)/JUN(N+1 )-ON/ALF 002190 
HN1=CMPLX(UN(N),-YN(N)) 002200 
RHON=HN1/HN-DN/ALF 002210 
AN®BN*S]GMAN/RHON 002220 
ANR(N)=REAL( AN) 002230 
ANI(N)#AIMAG(AN) 002240 

& R(N)=REAL( BN) 002250 

Bil (N)=ALTMAG( BN) 002260 
IF(IONETWO.EQ.2) GO TO 44 002270 
IF(NSCTOL.EQ.NTOL) GO TO 28 002280 
TERMSC=(2.0*DN+1.0)¢(ANR(N)*ANR(N)*ANI (N) SANI(N) 4 002290 

1 BNR(N)*BNR(N)+BNI(N)@BNI(N)) 002300 
CSCeCSC+TERMSC 002310 
NCSC=NCSC+1 002320 

iF (ABS( TERMSC/CSC).iT. TOL) 20,24 002330 

20 NSCTOL=NSCTOL+4 002349 
GO TO 28 002350 

24 NSCTOL=0 002360 
28 IF(NXTTOL.EQ-NTOL) GO TO 40 002370 
TERMXT=(2.0¢0N+1.0) *ANR(N)+8NR(N) 002380 
CXTECXT+TERMXT 002390 
NCXT=NCXT +1 002400 

IF (ABS( TERMXT/CXT). LT. TOL) 32,36 002410 

32 NATTOLENXTTOL +1 002420 
GO TO 40 002430 

36 NATTOL=0 002440 
40 CONTINUE 002450 
IF(NSCTOL.LT.NTOL.OR.NXTTOL.LT.NTOL) GO TO 44 002460 
NTAB=N 002470 

GO TO 48 002480 

44 CONTINUE 002490 
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NTAB*NHI 002500 

IF ( IONETWO.EQ.2) RETURN 002510 
CSCsCSC*LMDA*LMDA/(?.0*PI) 002520 
CXTsCXT*LMOASLMDA/(2.0¢PI1) 002530 
RETURN 002540 

END 002550 
002560 

002570 

002580 

SUBROUTINE SPHJN(X,N, ARR) 002590 
DIMENSION ARR(1) 002600 
RX21.£4270 002610 
JTS#AL0G10(RX)+0.5 002620 
OTS#ALOG10(RX) 002630 

NF=4 002640 
NN=N+1 002650 
1F(X.£Q.0.0) GO TO 4 002660 
IF(X.LE.0.5) GO TO 6 002670 

C FOR ORDERS LESS THAN ARGUMENT, AND ARGUMENT GREATER THAN 0.5 002680 
ARR(1)2#SIN(X)/X 002690 
ARR(2)*#ARR(1)/X-COS(X)/X 002700 
NM=X~2 002710 
1F(N.LT.X) NM=N-1t 002720 
1F(NM.LT.O) NM2O 002730 

DO 16 I=1,NM 002740 
FN=I+1 002750 

ARR( 142) =(FN+1.0)*ARR( I+?) /X-ARR(I) 002760 

C CHECK FO UNDERFLOW 002770 
IF(ARR(1+2).NE.0.0) GO TO 16 002780 
ARR(14+1) zARR( 141) RX 002790 
ARR(14+2)=(FN+1.0) #ARR( 141) /X-ARR( . ,4RX 002800 

16 CONTINUE 002810 
1F(N.LT.X) RETURN 002820 
HOLD=ARR(NM+2) 002830 

C FOR ORDERS GREATER THAN OR EQUAL TO ARGUMENT, AND 002840 
C. ARGUMENT GREATER THAN 0.5 002650 
15 ALPHA=1.0 002860 
M=~2-NM 002870 
NOTE=0 002880 
T=N+47 002890 
IF(X/N.GT.0.9) I#1+X/10.0 002900 
ARR(1+2)=ALPHA 002910 
ARR(1+3) 0.0 002920 
T=-I-1 002930 

DO 2 KK=1,M 002940 

K=-KK 002950 
NOO=K+4 002960 
TKPt=K+K+1 002970 
ARR(K)=TKP1*#ARR(K+1 ) /X-ARR(K+2) 002980 
OTS2=AL0G10({ TKP1 )+ALOG10(ASS(ARR(K+1) )) 002990 
IF(DTS2.GE.DTS.OR.OTS2-ALOG10(X).GE.DTS) 18,2 003000 

TEMPY =ARR(K) =ARR(K+1) /RX*¥ TKP1/X~ARR(K+2) /RX 003010 
ARR(K) =SARR(K) #RX 003020 
1F(ABS(ARR(K)).GE.1.0) 24,2 003030 
ARR(K)=TEMPY 003040 
ARR(K4+1) sARR(K+1) /RX 003050 
CONTINUE 003060 
C=HOLO/ARR(K) 003070 
IF(C.NE.0.0) GO TO 14 003080 
HOLD*HOLO*RX 003090 


NOTE=NOTE+1 003100 

GO TO 12 003110 

14 IF(NOTE.EQ.0) GG TO 17 003120 
MULT=JTS*NOTE 003130 

= 17 DO 3 I=K,NN 003140 
ARR(1)zARR(I)*C 003150 

3 CONTINUE 003160 
RETURN 003170 

C SPECIAL CASE WHERE X EQUALS 0.0 003180 
4 PRINT(6,%) 5 003190 
5 FORMAT(3X"*DO NOT USE ZERO ARGUMENT") 003200 
RETURN 003210 

C SERIES FOR ARGUMENT LESS THAN OR EQUAL To 0.5 003220 
6 EPS=SE~11 003230 
CH=CS=DN=DIV=1.0 003240 

00 10 L=1,NN 003250 
ITEXP=NEXPA=MEXPA=LEXPAaO 003260 
NO=L=4 003270 
IF(L.EQ.1) GO TO 7 003260 
CS=CS*X/DN 003290 

€ CHECK FOR UNDERFLOW 003300 
I1F(CS.NE.0.0) GO TO 7 003310 

C SCALE FOR THIS N AND ALL FOLLOWINF N 003320 
CS~CH*RX*¥X/DON 003330 

NE XPAsSNCXPAG1 003340 

7 DN=CN+2.0 003350 
CH=CS 003360 
XKK=0.0 003370 
A=CH2=FORTT21.0/DIV 003380 
TN=2.0eNO 003390 
FN=-0.5*x«X 003400 

8 XKK=XKK+1.0 003410 
FD=XKKe( TNtXKK+XKK44 0) 003420 
A=AeFN/FD 003430 

C CHECK FOR UNDERFLOW 003440 
1F(A.NE.0.0) GO TO 25 003450 

C SCALE FOR THIS N AND ALL FOLLOWING N 003460 
A=CH2*RX*FN/FD 003470 
O1V=DIV/RX 003480 
MEXPA=MEXPA+1 003490 
FORTT=FORIT#RX 003500 

2s Cri2=A 003510 
TT=A+FORTT 003520 

IF(TT .€Q. 0.0 .OR. ABS(A/TT) .LT. Eps) GO TO 9 003530 
FORTT=TT 003540 

GO TO 8 003550 

9 ARR(L)=CS*TT 003560 
C CHECK FOR UNDERFLOW 003570 
TFCARR(L).NE-0.0) GO TO 31 003580 

C SCALE FOR THIS N AND ALL FOLLOWING N 003590 
ARR(L)eCS#RX*TT 003600 
OlvsO1V/AXx 003610 
LEXPA=LEXPA+1 003620 

. 31 ITEXP=NEXPA+MEXPA+LEXPA 003630 
10 CONTINUE 003640 
END 003650 

c 003660 
c 003670 
Cc 003680 


SUBROUTINE SPHYN(X,N,ARR) 003690 


. 
a ote 


DIMENSION ARR(1) 003700 


Rx=1.€+270 003710 : 
DTS=ALOG?O( RX) 003720 er 
1F(X.EQ.0.0) 5,3 003730 aed 
PRINT(6,*) “OO NOT USE ZERO ARGUMENT FOR YN” 003740 ro 
RETURN 003750 “7 
FIND YO AND Y1 003760 << 
ARR(1)=-COS(X)/X 003770 ." 
ARR(2)=ARR(1)/X-SIN(X)/X 003780 <o 
NM1=N-1 003790 . * 
DO 2 K=1,NMt1 003800 
TKP1=K+K+1 003810 
ARR(K+2)=TKP1*ARR(K+t )/X-ARR(K) 003820 Be 
CHECK FoR UNDERFLOW 003830 eg 
OTS2=ALOGIO( TKPT )+AL0G10( ABS(ARR(K+1) )) 003840 ; 
IF (DTS2.GE.OTS.OR.DTS2-ALOG10(X).GE.OTS) 8,2 003850 Pixs 
ARR(K+1) fARR(K+t) /RX 003860 oes 
ARR(K+2)#TKP1*ARR(K+1)/X-ARR(K)/RX 003870 cae 
CONTINUE 003880 ., 
RETURN 003890 By 


END 003900 rs 
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ra Tables ae 
Zc 
This section contains a comprehensive set of theoretical scattering tables for 
2 = 
: conducting spheres. Each table lists normalized cross section(9/aa) for both aay 
‘ E-Plane and H-Plane at 2-degree intervals from backscatter (angle = 0°) to forward - 
scatter (angle = 180°) for a particular ka value. The ka values range from 0, 2 to ee 
: 40.0 at intervals of 0.2 and are accurate to at least four significant figures. ae 
- ~ 
aa 
a “. 
‘. “ 
ca Pa ie 
a 69 vu 
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MISSION 
of 
Rome Air Development Center 


RADC plans and executes research, development, test and 
selected acquisition programs in Support of Command, Control 
Communications and Intelligence (C71) activities. Technical 
and engineering support within aneas of technical competence 
44 provided to ESD Program Offices (POS) and other ESD 
elements. The principal technical mission areas are 
communications, electromagnetic guidance and control, sur- 
verllance of ground and aerospace ob jects, intelligence data 
collection and handling, Angormation system technology, 
sokid state sciences, electromagnetics and electronic 
reliability, maintainability and compatibivity. 
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